One-dimensional Bosons in Three-dimensional Traps 
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Bose gas ought to behave like the ID delta-function Bose gas solved by Lieb and Liniger. Up to 
now the theoretical arguments have been based on variational - perturbative ideas or numerical 
investigations. There are 4 parameters: density, transverse and longitudinal dimensions, and scat- 
tering length. In this paper we explicate 5 parameter regions in which various types of ID or 3D 
behavior occur in the ground state. Our treatment is based on a rigorous analysis of the many-body 
Schrodinger equation. 
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It appears to be possible to do experiments in highly 
elongated traps on ultra-cold Bose gases that are effec- 
tively ID. More precisely, the ID Bose gas with delta 
function 2-body interaction, analyzed long ago Q, 3, 
should be visible, with its quasi- fermionic behavior 3|, 
the absence of Bose-Einstein condensation (BEC) in a 
dilute limit 0. 0- @ ■ and an excitation spectrum differ- 
ent from that predicted by Bogoliubov's theory 0, • 
Several theoretical investigations on the transitions from 
3D to an effective ID behavior were triggered by bee, 
e.g., dmiUHElElHini. Systems showing indi- 
cations of such a transition have recently been prepared 
experimentally [II H EJ |H| . 

The theoretical work on the dimensional cross-over 
for the ground state in elongated traps has so far been 
based either on variational calculations, starting from a 
3D delta-potential I9lll5lll7l . or on numerical Quantum 
Monte Carlo studies |22L l23| with more realistic, genuine 
3D potentials, but particle numbers limited to the order 
of 100. This work is important and has led to valuable 
insi ghts , in particular about different parameter regions 
HoEll, but a more thorough theoretical understanding 
is clearly desirable since this is not a simple problem. In 
fact, it is evident that for a potential with a hard core the 
true 3D wave functions do not approximately factorize in 
the longitudinal and transverse variables (otherwise the 
energy would be infinite) and the effective ID potential 
can not be obtained by simply integrating out the trans- 
verse variables of the 3D potential (that would immedi- 
ately create an impenetrable barrier in ID). It is impor- 
tant to be able to demonstrate rigorously, and therefore 
unambiguously, that the ID behavior really follows from 
the fundamental Schrodinger equation. It is also impor- 
tant to delineate, as we do here, precisely what can be 
seen in the different parameter regions. The full proofs of 
our assertions are long and will be given elsewhere |24j . 



but we emphasize that everything can be rigorously de- 
rived from first principles. In this paper we state our 
main results and outline the basic ideas for the proofs. 

We shall always be concerned with the ground state 
and with large particle number, N 1, which is appro- 
priate for the consideration of actual experiments. In or- 
der to make precise statements we take the limit N — > oo 
but the reader can confidently apply these limiting state- 
ment to finite numbers like N = 100. Other parameters 
of the problem are the scattering length, a, of the two- 
body interaction potential, v, and two lengths, r and L, 
describing the transverse and the longitudinal extension 
of the trap potential, respectively. 

It is convenient to write the Hamiltonian in the follow- 
ing way (in units where h — 2m = 1): 

N 
3 = 1 

+ WaGxi-X^I) (1) 



l<i<j<N 



with x = (x, y, z) = (x^, z) and with 



Vl(z) 



L 2 



K x (x x ) = ^y ± (x», 

V(z/L), Wo (|x|) = l«(|x|/ a ) . (2) 



Here, r,L,a are variable scaling parameters while V , 
V and v are fixed. The interaction potential v is sup- 
posed to be nonnegative, of finite range and have scatter- 
ing length 1; the scaled potential v a then has scattering 
length a. The external trap potentials V and V confine 
the motion in the longitudinal (z) and the transversal 
(x- 1 ) directions, respectively, and are assumed to be con- 
tinuous and tend to oo as \z\ and |x^| tend to oo. To 
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simplify the discussion we find it also convenient to as- 
sume that V is homogeneous of some order s > 0, namely 
V(z) = | z I s , but weaker assumptions, e.g. asymptotic ho- 
mogeneity |25j, would in fact suffice. The case of a simple 
box with hard walls is realized by taking s = oo, while 
the usual harmonic approximation is s = 2. It is under- 
stood that the lengths associated with the ground states 
of -d 2 /dz 2 + V(z) and -(V ± ) 2 + V ± (x L ) are both of the 
order 1 so that L and r measure, respectively, the lon- 
gitudinal and the transverse extensions of the trap. We 
denote the ground state energy of by E^ M (N, L,r,a) 
and the ground state particle density by p ( ff^ r a (x). 

While the 3D density is always assumed to be low (in 
the sense that distance between particles is large com- 
pared to the 3D scattering length) the ID density can be 
either high or low. In contrast to 3D gases, high density 
in ID corresponds to weak interactions and vice versa 0. 

To keep the discussion simple let us first think of the 
case that the particles are confined in a box with di- 
mensions r and L. The 3D particle density is then 
p 3D = N/{r 2 L) and the ID density p 1D = N/L. The 
case of quadratic or more general homogeneous trapping 
potentials will be considered later. We begin by describ- 
ing the division of the space of parameters into two basic 
regions. This decomposition will eventually be refined 
into five regions, but for the moment let us concentrate 
on the basic dichotomy. 

In earlier work [2jJ |26j we proved that the 3D Gross- 
Pitaevskii formula for the energy (including its limiting 
'Thomas-Fermi' case) is correct to leading order in situ- 
ations in which a is small and N is large. This energy 
has two parts: The energy necessary to confine the parti- 
cles in the trap, which is roughly N(r~ 2 +L~ 2 ), plus the 
internal energy of interaction, which is iV47rap 3D . The 
trouble is that while this formula is correct for a fixed 
confining potential in the limit N — > oo with a 3 p 3D — » 0, 
it does not hold uniformly if r/L gets small as N gets 
large. In other words, new physics can come into play as 
r/L — ► and it turns out that this depends on the ratio 
of a/r 2 to p 1D — N/L . There are two basic regimes 
to consider in highly elongated traps, i.e., when r <C L. 
They are are 

• The ID limit of the 3D Gross-Pitaevskii/'Thomas- 
Fermi' regime 

• The 'true' ID regime. 

The former turns out to be characterized by aL/r 2 N — > 
0, while in the latter regime aL/r 2 N is of the order one 
or even tends to infinity. (This is usually referred to as 
the Girardeau- Tonks region, but it was Girardeau 
not Tonks, who understood how to calculate the states 
of the ID hard-core gas.) These two situations corre- 
spond to high ID density (weak interaction) and low ID 
density (strong interaction), respectively. Physically, the 
main difference is that in the strong interaction regime 



the motion of the particles in the longitudinal direction 
is highly correlated, while in the weak interaction regime 
it is not. Mathematically, this distinction also shows up 
in our proofs. 

In both regions the internal energy of the gas is small 
compared to the energy of confinement which is of the 
order N/r 2 . However, this in itself does not imply a 
specifically ID behavior. (If a is sufficiently small it is 
satisfied in a trap of any shape.) ID behavior, when it 
occurs, manifests itself by the fact that the transverse 
motion of the atoms is uncorrelated while the longitudi- 
nal motion is correlated (very roughly speaking) in the 
same way as pearls on a necklace. Thus, the true crite- 
rion for ID behavior is that aL/r 2 N is of the order unity 
or larger, and not merely the condition that the energy 
of confinement dominates the internal energy. 

In parallel with the 3D Hamiltonian we consider the 
Hamiltonian for n Bosons in ID with delta interaction 
and coupling constant g > , i.e., 



We consider this Hamiltonian for the Zj in an interval of 
length £ in the thermodynamic limit, £ — ► oo, n — ► oo 
with p = n/£ fixed. The ground state energy per particle 
in this limit is independent of boundary conditions and 
can, according to |l| , be written as 



e^(p)=p J e(g/p) , 



(4) 



with a function e(t) determined by a certain integral 

3 

and e(t) — » 7r 2 /3 for t — > oo. Thus 



equation. Its asymptotic form is e{t) as ^t for { < 1 



-ID 



(p) ~ hp for g/p < 1 



and 



el B (p) « (tt7% 2 for 5 /p»l 



(5) 



(6) 



Taking pe}p{p) as a local energy density for an inho- 
mogeneous ID system we can form the energy functional 



£\P\ = / (|V^(z)| 2 + V L (z)p(z) + p(zfe{g/p{z))) dz 



with ground state energy defined to be 
E 1D (N, L,g) = inf <£[p] : p(z)>0, [ p(z)dz = N 



(7) 



(8) 

By standard methods (cf., e.g., 26]) one can show that 
there is a unique minimizer, i.e., a density PN,L,g(z) with 

/ PN, L A z ) dz = N and £\pN,L,g] = E 1D (n\L i9 ). We 

define the mean ID density of this minimizer to be 



P=Jf I {PN,L,g{ Z )f dz 



(9) 
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In a rigid box, i.e., for s = oo, p is simply N/L (except for 
boundary corrections), but in more general traps it de- 
pends also on g besides N and L. The order of magnitude 
of p in the various parameter regions will be described 
below. 

Our main result relates the 3D ground state energy of 
JTJ, E QM (N, L, r, a), to the ID density functional energy 
E 1D (N, L,g) in the large N limit with g ~ a/r 2 provided 
r/L and a/r are sufficiently small. To state this precisely, 
let e and ^(x- 1 -), respectively, denote the ground state 
energy and the normalized ground state wave function 
of — (V^) 2 + y ± (x ± ). The corresponding quantities for 
-(V ± ) 2 +V;- L (x ± ) aree-L/r 2 and6 r (x- L ) = (l/r)6(x ± /r). 
In the case that the trap is a cylinder with hard walls b 
is a Bessel function; for a quadratic V 1 - it is a Gaussian. 

Define g by 



|6(x x )| 4 d 2 x x =8na / |6 r (x- L )| 4 d 2 x- L . (10) 



Theorem. Let N — > oo and simultaneously r/L — > 
and a/r — > in such a way that r 2 p ■ min{p, g} — ► 0. 
TTien 



lim -fr- : = 1. 

E™(N,L,g) 



(11) 



Moreover, if we define the ID quantum mechanical den- 
sity by averaging over the transverse variables, i.e., 



PN,L,r,a( Z ) = 



(12) 



then p]^ r a ( z )/ Pn,l,q(z) — »• 1 in a suitable sense. 



Note that because of (JSJ) and © the condition r 2 p 
min{p, g} — > is the same as 



4 D (p) « 1/r 2 , 



(13) 



i.e., the average energy per particle associated with the 
longitudinal motion should be much smaller than the en- 
ergy gap between the ground and first excited state of the 
confining Hamiltonian in the transverse directions. Thus, 
the basic physics is highly quantum-mechanical and has 
no classical counterpart. The system can be described 
by a ID functional Q, even though the transverse trap 
dimension is much larger than the range of the atomic 
forces. 

The two regimes mentioned previously correspond to 
specific restrictions on the size of the ratio g/p as N — > 
oo, namely g/p <C 1 for the limit of the 3D Gross- 
Pitaevskii regime (weak interaction/high density), and 
g/p > for the 'true' ID regime (strong interaction/low 
density). The precise meaning of <C is that the ratio of 
the left side to the right side tends to zero in the limit 
considered. 

We shall now briefly describe the finer division of these 
two regimes into five regions altogether. Three of them 



(Regions 1-3) belong to the weak interaction regime and 
two (Regions 4-5) to the strong interaction regime. In 
each of these regions the general functional Q can be 
replaced by a different, simpler functional, and the energy 
E 1B (N,L,g) in the theorem by the ground state energy 
of that functional, in analogy with the fact that in 3D the 
TF functional, i.e., the GP functional without gradient 
term, is a limit of the full GP functional. We state the 
results for the ground state energy, but corresponding- 
results hold for the density. We always assume N — ► oo 
and r/L — > 0. 

The five regions are 

• Region 1, the Ideal Gas case: g/p <C ^V~ 2 , with 
p r*> N/L, corresponding to the trivial case where the 
interaction is so weak that it effectively vanishes in the 
large N limit and everything collapses to the ground state 
of —d 2 /dz 2 + V{z) with ground state energy e". The 
energy E 1D in l|ll(l can be replaced by Ne^/L 2 . Note 
that g/p iV~ 2 means that the 3D interaction energy 
- j0 3]i a« l/L 2 . 

• Region 2, the ID GP case: g/p ~ iV -2 , with 
p ~ N/L, described by a ID Gross-Pitaevskii energy 
functional 

/oo 
(\\7^(z)\ 2 + V L (z)p(z) + y P (z) 2 ) dz 
-OO 

(14) 

corresponding to the high density approximation (J5J of 
the interaction energy in 0. Its ground state energy is 
E GP (N,L,g) = NL- 2 E GP (1,1, NgL), by scaling. 

• Region 3, the ID TF case: N~ 2 < g/p < 1, with 
p ~ (N/L)(NgL)^ 1 ^ s+1 \ where s is the degree of homo- 
geneity of the longitudinal confining potential V. This 
region is described by a Thomas-Fermi type functional 



£ TF [P}= (V L {z)p{z) + \gp{z) 2 )dz 



(15) 



It is a limiting case of Region 2 in the sense that 
NgL ~ NaL/r 2 — * oo, but a/r is sufficiently small 
so that g/p - (aL/ATr 2 )(ATaL/r 2 ) 1 /( s + 1 ) -> 0, i.e., 
the high density approximation in JSJl is still valid. In 
this limit the gradient term in l|14|) becomes vanish- 
ingly small compared to the other terms. The expla- 
nation of the factor {NaL/r 2 ) 1 /^ 1 ^ - {NgL) 1 H s+l ) 
is as follows: The linear extension L of the minimiz- 
ing density p^v 1 ^ g is for large values of NgL determined 
by V L {L) - g(N/L), which gives L - {NgLf/^+^L. 
In addition condition l|13(l requires gp r~ 2 , which 
means that Na/L{NaL/r 2 ) 1 / ( - s+1 ^ -> 0. The minimum 
energy of l|15|l has the scaling property E TF (N,L,g) = 
NL- 2 {NgL) s ^ s +^E TF (1,1,1). 

• Region 4, the LL case: g/p ~ 1 , with p ~ 
{N/L)N- 2 / {s+2 \ described by an energy functional 



(V L (z)p(z)+p(z) 3 e(g/p(z)))dz . (16) 
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This region corresponds to the case g/p ~ 1, so that 
neither the high density JSJ nor the low density approxi- 
mation (JSJ is valid and the full LL energy Q has to be 
used. The extension Z of the system is now determined 
by V L {L) ~ (N/L) 2 which leads to Z LL = LN 2 ^ s+2 \ 
in contrast to L TF = L(NgL) 1 / ( - s+r > in the TF case. 
Condition l|13|) means in this region that iW/Zjx ~ 
N s/(s+2) r i L _^ Q Since 7Vr/Z LL - (p/g)(a/r), this con- 
dition is automatically fulfilled if g/p is bounded away 
from zero and a/r —>■ 0. The ground state energy of (|16J) . 
£ LL (N, L,g), is equal to iV7 2 £: LL (l, 1, g/-y), where we 
introduced the density parameter 

7 = A/Z LL = {N/L)N- 2 l {s+2 *> . (17) 

• Region 5, the GT case: g/p ^ 1, with p ~ 
(N/L)N~ 2 /( s+2 \ described by a functional with energy 
density ~ p 3 , corresponding to the Girardaeu-Tonks limit 
of the the LL energy density. It corresponds to impene- 
trable particles, i.e, the limiting case g/p — > oo and hence 
formula JJJJ for the energy density. As in Region 4, the 
mean density is here p ~ 7. The energy functional is 

£ GT M = I" (Vl(z)p(z) + (n 2 /3)p(zf) dz , (18) 
J — 00 

with minimum energy E GT (N, L) — Nj 2 E GT (l, 1). 

We note that the condition g/p ~ 1 means that Region 
4 requires the gas cloud to have aspect ratio r/L of the 
order N^ 1 (a/r) or smaller, where Z ~ LN 2 /( S + 2 ) is the 
length of the cloud. Experimentally, such small aspect 
ratios are quite a challenge and the situations described 
in |l8lll9l l20ll21| are still rather far from this regime. It 
may not be completely out of reach, however. 

Regions 1-3 can be reached as limiting cases of a 3D 
Gross-Pitaevskii theory. In this sense, the behavior in 
these regions contains remnants of the 3D theory, which 
also shows up in the the fact that Bose-Einstein conden- 
sation (BEC) prevails in Regions 1 and 2 j^. Heuris- 
tically, these traces of 3D can be understood from the 
fact that in Regions 1-3 the ID formula for energy per 
particle, gp ~ aN/(r 2 L), gives the same result as the 
3D formula jsjj, i.e., scattering length times 3D density. 
This is no longer so in Regions 4 and 5 and different 
methods are required. 

BEC probably occurs in part of Region 3, besides Re- 
gions 1 and 2, but we cannot prove this and it remains an 
open problem. BEC means [2JJ that the one-body den- 
sity matrix 7(x, x') factorizes as iV</>(x)0(x') for some 
normalized (j>. This, in fact, is 100% condensation and 
this is what we prove occurs (in the N — ► 00 limit, of 
course). The function (f> is the square- root of the mini- 
mizer of the ID GP functional (|14|l times the transverse 
function ^(x- 1 ). The proof is similar to the work in (271 • 
BEC is not expected in Regions 4 and 5. Lenard |4| 
showed that the largest eigenvalue of 7 grows only as 



iV 1 / 2 for a homogeneous gas of ID impenetrable bosons 
and, according to [28| and [2^, this holds also for a GT 
gas in a harmonic trap. (The exponent 0.59 in can 
probably be ascribed to the small number of particles 
(N = 10) considered.) 

We now outline the main steps in the proof of the The- 
orem, referring to [24| for full details. The different pa- 
rameter regions have to be treated by different methods, 
a watershed lying between Regions 1-3 on the one hand 
and Regions 4-5 on the other. In Regions 1-3, similar 
methods as in the proof of the 3D Gross-Pitaevskii limit 
theorem in |2g can be used. This 3D proof needs consid- 
erable modifications, however, because in (2(| the exter- 
nal potential is fixed and the estimates are not uniform 
in the ratio r/L. 

To prove (|ll|l one has to establish upper and lower 
bounds, with controlled errors, on the QM many-body 
energy in terms of the energies obtained by minimiz- 
ing the energy functionals appropriate for the various re- 
gions. The limit theorem for the densities can be derived 
from the energy estimates in a standard way by varia- 
tion with respect to the external potential Vl ■ As usual, 
the upper bounds for the energy are easier than the lower 
bounds, but nevertheless not simple, in particular not for 
'hard' potentials v. 

The upper bound in Regions 1-3 is ob- 
tained from a variational ansatz of the form 
*(xi,...,Xjv) = F(Ki,...,x N )H^ =1 b r (x^)y/p^ w (z k ), 

with F(xi, . . . , Xjv) = rifc=i /( x fc - x j(fc)) where x j(fc) is 
the nearest neighbor of x^ among the points Xj, j < k, 
b r (x- L ) is the lowest eigenfunction of — (V^) 2 + V r , and 
p GP (z) the minimizer of the ID GP functional 1)14(1 . The 
function / is, up to a cut-off length that has to be chosen 
optimally, the zero energy scattering solution for the 
two-body Hamiltonian with interaction v a . The form 
of F, inspired by [3]]], is chosen rather than a Jastrow 
ansatz Yii<j /( x « ~ x j) because it is computationally 
simpler for the purpose of obtaining rigorous estimates. 

For an upper bound in Regions 4-5 a natural vari- 
ational ansatz would appear to be \&(xi, . . . , xjv) = 

F(xi, . . . ,Xjv)rifcLA( x fc ■ ■ ■ ,zn) where ip is the 

ground state of Hjf g with the external potential Vl 
added. However, in order to make a link with the exact 
solution Q for a homogeneous gas, but also to control 
the norm of the trial function, it turns out to be nec- 
essary to localize the particles by dividing the trap into 
finite 'boxes' (finite in z-direction), with a finite parti- 
cle number in each box and making the ansatz with the 
boundary condition ^> = for each box individually. The 
particles are then distributed optimally among the boxes 
to minimize the energy. This box method, but with the 
boundary condition = 0, is also used for the lower 
bounds to the energy. Another essential device for the 
lower bounds is Dyson's Lemma that was also used in 
HU E3- H3- This lemma, which goes back to Dyson's 



5 



seminal paper [31( on the hard core Bose gas, estimates 
the kinetic and potential energy for a Hamiltonian with 
a 'hard' potential v of finite range from below by the po- 
tential energy of a 'soft' potential U of larger range but 
essentially the same scattering length as v. Borrowing a 
tiny part of the kinetic energy it is then possible to do 
perturbation theory with the 'soft' potential U and use 
Temple's inequality to bound the errors. A direct ap- 
plication of perturbation theory to the original potential 
v, on the other hand, is in general not possible. 

A core lemma for Regions 4-5 is a lower bound on the 
3D ground state energy in a finite box in terms of the ID 
energy of the Hamiltonian both with the boundary 
condition V\& = 0. Denoting the former energy by E^ 
and the latter by E^ K , this bound for n particles in a 
box of length £ in the z-direction reads 



^box ' 



1 - Cn 



nr f a\ V 8 



with a constant C. To prove this bound the ground state 
wave function is first written as a product of JJ. b r (xj:) 
and a function G(xi, . . . , x„). This subtracts ne^/r 2 
from E^ but the resulting minimization problem for G 
involves the weighted measure JJ k b r (x^) 2 d 3 Xk in place 
of rife^ 3x fc- Nevertheless, Dyson's Lemma can be used, 
and after the 'hard' potential v has been replaced by the 
'soft' potential U, it is possible to integrate the trans- 
verse variables away and obtain a minimization prob- 
lem for a ID many body Hamiltonian with interaction 
d{ Zl - Zj ) = a J ^6 r (x^) 2 6 r (x|) 2 t/(x i - Xj-JdVd 2 ^. 
In the limit considered this converges to a delta interac- 
tion with the coupling constant IjlOfl. The error terms in 
the estimate for E^ x arise both from Temple's inequal- 
ity and the replacement of d by a delta function, among 
other things. When the particles are distributed opti- 
mally among the boxes to obtain a global lower bound, 
superadditivity of the energy and convexity of the energy 
density p 3 e(a/ p) are used, generalizing corresponding ar- 
guments in |2flj. 

In conclusion, we have reported a rigorous analysis of 
the parameter regions in which a Bose gas in an elongated 
trap may or may not be expected to display ID behav- 
ior in its ground state. This takes the form of theorems 
about the ground state energy and density. We always 
consider N to be large and r/L to be small, but this is 
not enough to distinguish ID from 3D behavior. The dis- 
tinction occurs when, additionally, the ID density is not 
too large, e.g., aL/r 2 ~ N 3 / 2 in a harmonic trap. The 
ID behavior, when it occurs, is described by the delta- 
function gas solved in [jj. We also present a ID energy 
functional, analogous to the Gross-Pitaevskii functional, 
that correctly describes the energy and density in all the 
5 parameter regions considered here. 
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